Initial state dependence of the quench dynamics in integrable quantum systems 
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We identify and study classes of initial states in integrable quantum systems that, after the 
relaxation dynamics following a sudden quench, lead to near-thermal expectation values of few- 
body observables. In the systems considered here, those states are found to be insulating ground 
states of lattice hard-core boson Hamiltonians. We show that, as a suitable parameter in the 
initial Hamiltonian is changed, those states become closer to Fock states (products of single site 
states) as the outcome of the relaxation dynamics becomes closer to the thermal prediction. At 
the same time, the energy density approaches a Gaussian. Furthermore, the entropy associated 
with the generalized canonical and generalized grand-canonical ensembles, introduced to describe 
observables in integrable systems after relaxation, approaches that of the conventional canonical 
and grand-canonical ensembles. We argue that those classes of initial states are special because a 
control parameter allows one to tune the distribution of conserved quantities to approach the one in 
thermal equilibrium. This helps in understanding the approach of all the quantities studied to their 
thermal expectation values. However, a finite-size scaling analysis shows that this behavior should 
not be confused with thermalization as understood for nonintegrable systems. 
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I. INTRODUCTION 

The relaxation dynamics of isolated quantum systems 
after a sudden quench is a topic that is attracting much 
current attention. Interest on this problem has been 

tarkcd by recent experiments with ultracold gases [l|- 
The high degree of isolation in those experiments 
allows one to consider them as almost ideal analog sim- 
ulators of the unitary dynamics of pure quantum states. 
For example, in Ref. [4I, Kinoshita et al. showed that 
observables in a (quasi-)one-dimensional bosonic system 
close to an integrable point do not relax to the values 
expected from a conventional statistical mechanics de- 
scription. Any non-negligible coupling to a thermal en- 
vironment would have destroyed such a remarkable phe- 
nomenon. More recently, Trotzky et al. [J| have shown 
that the experimental dynamics of Bose-Hubbard like 
(quasi-)one-dimensional systems can be almost perfectly 
described by the unitary dynamics of the relevant model 
Hamiltonian. The latter was followed by numerically ex- 
act means utilizing the time-dependent renormalization 
group algorithm 0, @ . 

After the experimental results in Ref. [4I, many theo- 
retical works have found that, following a sudden quench 
within integrable systems, few-body observables, in gen- 
eral, relax to nonthermal expectation values J7H23| (for a 



recent review, see Ref. 12J]). Some of the novel insights 
gained through these studies include (i) the possibility of 
describing observables after relaxation by means of gen- 
eralized Gibb ensembles (GGE) 0-[ll [mH EHll ; (ii) 
the fact that even though in some cases the behavior of 
nonlocal observables after relaxation can be parametrized 
similarly to the one in thermal equilibrium |15l . |20[ , an 
exact description of those observables is provided only 
by the GGE Q, and (iii) an understanding of the GGE 
through a generalization of the eigenstate thermalization 



hypothesis (ETH) [21|. ETH explains why thermaliza- 
tion occurs in generic (nonintegrable) quantum systems 
after a quench |25l - [27| . 

All the results discussed above have been obtained in 
studies of several specific models. However, they are ex- 
pected to hold in general for integrable systems. An 
interesting, and so far nongeneric, result reported in 
Refs. [I,l21| was the observation of a phenomenon close to 
"real" thermalization in integrable systems, in the sense 
of the expectation values of few-body observables after 
relaxation approaching those predicted in thermal equi- 
librium. This occurred as a parameter used to generate 
special classes of initial states was changed. In Ref. Q , 
the initial states were insulating ground states of hard- 
core bosons in half-filled period-two superlattices, while 
in Ref. [21| , they were the ground state of trapped sys- 
tems with a Mott insulating domain in the trap center. 

In this work, we revisit the systems above and fo- 
cus on understanding the properties of the initial states 
for which observables after relaxation were seen to ap- 
proach thermal expectation values, despite integrability. 
As said before, those states are insulating ground states. 
Here, we show that the selected tuning parameter makes 
those initial states approach Fock states (products of sin- 
gle site wavefunctions) at the same time that (i) their 
energy density approaches a Gaussian, and (ii) the en- 
tropy of their associated generalized canonical and grand- 
canonical ensembles approach the entropies of the canoni- 
cal and grand-canonical ensembles. We argue that (i) and 
(ii) above can be understood because the distribution of 
conserved quantities in such initial states approaches the 
one of systems in thermal equilibrium. Hence, they can 
have thermal-like energy densities, entropies, and observ- 
ables after relaxation. However, after a finite-size scal- 
ing analysis, we conclude that this phenomenon differs 
conceptually from thermalization as it happens in noin- 



tegrable systems. 

The presentation is organized as follows: in Sec. HIl 
we introduce the models and observables of interest. We 
also define the ensembles considered and provide details 
on how the calculations are performed. Section |lll] is de- 
voted to study of the overlaps of the initial states with 
the eigenstates of the final Hamiltonians, as well as to 
the description of the energy densities in all cases. The 
scaling of the entropy with system size, for the different 
ensembles analyzed and for superlattice and trapped sys- 
tems, is presented in Sec. HV] In Sec. El we study the 
distribution of the conserved quantities for the different 
initial states and within standard statistical ensembles. 
Finally, the conclusions are presented in Sec. |Vll 



II. MODEL, ENSEMBLES, AND OBSERVABLES 

We are interested in the equilibrium and noncquilib- 
riuni properties of lattice bosons in the limit of infinite 
on-site repulsion (hard-core bosons). Those systems can 
be described by the Hamiltonian 



the dynamics, which satisfies ifj^'a) = Ea\"^a)- Hence, 
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where D is the dimension of the Hilbert space and 
Ca = {^a\'>Pi)^ s-nd the time evolving wave function can 
be written as 



D 



|*(i)) = e-'"'/^\^i) = J2 ae-^^°*/''|vI'„). (3) 
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The time evolution of a generic observable O is then 
dictated by the sums over all eigenstates 

(0(t)) = (vl'(t)IOIvI'(i)) =^c:C0e^(^"-^^)*/''O„^, 

a,l3 

where Oap = (^a|0|^/3), and the infinite time average 
of Eq. (m (in the absence of degeneracies) can be thought 
as the result of a diagonal ensemble average [23| 



(0)de = ^|C„|20, 
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As shown in Refs. [21|, |27[, this infinite time average de- 
scribes observables after relaxation. We should stress 
that this can be true even in the presence of degenera- 
cies associated with i nteg rability, except for cases with 



with the additional on-site constraints ir- = M = 0, 

which preclude multiple occupancy of the lattice sites. massive degeneracies [Ij. The" validity of the descrip- 
Here, J is the nearest-neighbor hopping, V-"^^ is a site- 
dependent local potential, and L is the number of lattice 
sites. The hard-core boson creation (annihilation) opera- 



tor in each site is denoted by 6 J {hj) and the site number 

occupation by fij = ijbj- In what follows, we consider 
only systems with open boundary conditions, and f = 1 
sets our units of energy. 

This model is integrable [25 and can be exactly solved 
by first mapping it onto the spin- 1/2 XX model (with 
a site-dependent magnetic field in the z direction) by 
means of the Holstein-Primakoff transformation [2^ and 
then onto a noninteracting fcrmion Hamiltonian utiliz- 
ing the Jordan- Wigner transformation [28, 30i | . In the 
fermionic language, the Hamiltonian can be straightfor- 
wardly diagonalized and all the spectral and thermody- 
namic properties of hard-core bosons can be computed 
either analytically or numerically in polynomial time. 
Off-diagonal correlations are more difficult to calculate. 
However, using properties of Slater determinants, they 
can also be computed very efficiently numerically for 
ground state [3l|, [S^ and finite temperature [SjI equi- 
librium problems, as well as during the unitary nonequi- 
librium dynamics j34| . Those insights will be used later. 

More generally, the nonequilibrium dynamics of iso- 
lated quantum systems can be studied by writing the 
(arbitrary) initial state \^i) as a linear combination of 
the eigenstates J^q) of the Hamiltonian H that drives 



and for the same 
23 



tion of integrable systems after relaxation, by means of 
the infinite time average ([5]), has been demonstrated for 
the l/r Hubbard model in Ref. [^ 
(hard-core boson) systems considered here in Ref. 
(supplementary materials) . 

The result in Eq. ([S]) is to be compared with the pre- 
dictions of conventional statistical mechanics ensembles, 
for a system in equilibrium with energy Ei = {ipi\H\il)i) 
and total number of particles N . The canonical ensemble 
predicts 
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where Zqe — X^a ^ Ea/ksT ^ rp j^gg^jg ^q j^g ^^ken such 
that Ej = Z^^ J2a e~^°/'^^"^i?Q, and the sums run over 
all eigenstates of the Hamiltonian (with energy Ea) in the 
sector with N particles. The grand-canonical ensemble, 
on the other hand, predicts 
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GE 
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A^' 






(7) 



where Zqe ^ J2o 



-(£„ 



-f^N^)/kBT^ T and ^ need to be 



taken such that Ei = Z^^J2o.e^^^''~^^"''^''''^ E^ and 
^ = ^GE Ea e"(^°-''^<»)/'=^'^7V„, and the sums run over 
all eigenstates of the Hamiltonian (with energy Ea and 



number of particles Na). The predictions of Eq. ((6)) and 
Eq. ([7]) in general agree in the thermodynamic limit. 

Thermalization is then said to occur if, for sufficiently 
large systems, (0)de — (0)ce — (0)ge- Hence, the fact 
that thermalization occurs in isolated systems is surpris- 
ing as (0)de depends on the initial conditions through 
the projection of the initial state onto all the eigenstates 
of the Hamiltonian, while conventional statistical ensem- 
bles depend only on the initial conditions through Ej 
and N. Since the energy distribution of the initial state 
li/jj) in the eigenstates of the final Hamiltonian is nar- 
row (because of locality, see Ref. [231 ^^^ its supplemen- 
tary materials) and centered around Ej , similarly to the 
canonical and grand-canonical ensembles, then thermal- 
ization can be understood to occur because of ETH |25l - 
[271 . ETH states that in generic many-body systems Oaa 
almost do not fluctuate between eigenstates that have 
similar energies, i.e., the eigenstates themselves already 
exhibit thermal behavior. 

Furthermore, it has been also proposed that one can 
define the entropy of the isolated system after the quench 
to be the diagonal entropy j35| 



Sd^-J2\Cc.\''H\Co.\ 



(8) 



which satisfies all the thermodynamic properties required 
from an entropy. Indeed, this entropy has been recently 
shown to be consistent with the microcanonical entropy 
for nonintegrable systems [3g, and hence, for sufficiently 
large systems it is expected to agree with the entropy of 
the canonical ensemble 
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and with that of the grand-canonical ensemble 
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(10) 



up to subextensive corrections. 

In general, in integrable systems such as the ones of in- 
terest in this work, the presence of a complete set of con- 
served quantities prevents thermalization [3, [231 ■ (The 
eigenstate thermalization hypothesis has been shown to 
fail in those systems j2ll.[27|.) However, after relaxation, 
few-body observables can be described by means of a 
generalization of the Gibbs ensemble |7[ , with a density 
matrix 



Pgge 



Z, 



where ^gge = Tr 



GGE^ 



-E„A„/„ 



-E„A„/„ 



(11) 



{/„} are the conserved 



quantities, and n = 1, . . . , L. In our systems, {/«} are the 
occupation operators of the single-particle eigenstates of 
the noninteracting fermionic Hamiltonian to which hard- 
core bosons can be mapped. {A„} are the Lagrange 
multipliers, which are selected such that {ipj\In\ipi) — 



Tr(J„pGGE). For hard-core bosons, they can be com- 
puted using the expression [31 



(12) 



A„ = In 


' I - {lPl\in\^Pl)' 
(V'7|4|V'7) 


^GGE is then 


^GGE = 


=n(i+e-'")- 



(13) 



The fact that the GGE is able to predict expectation 
values of few-body observables after relaxation can be un- 
derstood in terms of a generalized ETH [2l| . The idea in 
this case is that eigenstates of the Hamiltonian that have 
similar values of the conserved quantities have similar 
expectation values of few-body observables. The GGE 
is then the ensemble that, within the full spectrum, se- 
lects a narrow set of states with the same distribution of 
conserved quantities that is fixed by the initial state. 

The GGE entropy is given by 



•^GGE = IuZgGE + ^ Xn{lpl\In\lpl)- 



(14) 



Furthermore, one can also define a canonical version of 
this generalized ensemble, with a density matrix 



-7-1 -V ^-ni-n 

PGCE - ^GCE^' " ' 



(15) 



for which only states with N particles are considered 
when calculating traces. We keep A„ in the sector with N 
particles to have the same values as within the GGE and 
take the partition function to be the trace over states 

The en- 



with N particles, Zqce = Tr e"^.. ^"^" 
tropy of this ensemble can be computed as 



N 



•S'gce = Tr[/5GCE In(pGCE)] 



N 



(16) 



where, once again, only eigenstates of the Hamiltonian 
with N particles contribute to the trace. 

An interesting recent finding in Ref. [36[ was that de- 
spite the fact that the generalized ensembles do describe 
few-body observables in integrable systems after relax- 
ation, their entropy is always greater than that of the 
diagonal ensemble, and the difference increases linearly 
with increasing system size. This means that an expo- 
nentially larger number of states contribute to the gen- 
eralized ensembles when compared to the diagonal one. 
The generalized ETH ensures that, despite having a much 
greater number of states, the generalized ensembles pre- 
dict the outcome of the realization dynamics. This is 
because the overwhelming majority of the states that 
contribute to the generalized ensembles have identical 
expectation values of few-body observables as the ones 
that contribute to the diagonal ensemble [2l|. All these 
results are expected to be generic in integrable systems. 

In this work, instead, we focus on special classes of 
initial states that lead to expectation values of few-body 
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FIG. 1: (Color online) Weights of the eigenstates of the final 
Hamiltonian in the diagonal (top half in both panels) and 
canonical (bottom half in both panels) ensembles, |Cap and 
e"^°/*«^/ZcE, respectively, for L = 36 a.nd N = 18 (half 
filling). The panel on the left (a) depicts results for a quench 
from Ai = 4 to Af ~ 0, and the panel on the right (b) for a 
quench from Ai = to Af — 4. In both cases, we select the 
initial state to be the ground state of Eq. ^ for the given 
value of A = Aj. The color scale indicates the number of 
states, per unit area in the plot, that have a given weight. 



observables that approach those in thermal equilibrium, 
despite integrability [^, |2l| . Since we know that ETH is 
not satisfied in those systems [2l|, [23| , the fact that ob- 
servables after relaxation approach thermal values then 
must be related to special properties of the overlaps Cq of 
the initial states with the eigenstates of the final Hamil- 
tonians. Hence, we study the behavior of the Cq,'s in such 
systems. Hard-core bosons can be mapped onto nonin- 
teracting fermions, so one can generate the exponentially 
large Hilbcrt space of finite systems [whose size is (^)] 
without the need of diagonalizing the full Hamiltonians. 
Those many-body states are created as products of non- 
interacting fcrmionic eigenstates. They can be written as 
Slater determinants, in terms of fcrmionic creation opcr- 
t 



ators fl, as 



N L 

1=1 fc=i 



(17) 



1^/) - 



and the same can be done for the initial state 

Yli=iJ2k=i^ki fk\^)- "^^^ overlap between the initial 
state and the eigenstates of the final Hamiltonian can 
then be calculated numerically as the determinant of the 
product of two matrices [3l|, |33| 



C„ = (*o|^/)=det (P")^P" 



(18) 



FIG. 2: (Color online) Same as Fig. [T] but for L = 44 and 
Af = 11 (quarter filling). 



which, together will all the expressions presented pre- 
viously, allow us to compute the energy distributions 
and entropies in the diagonal, canonical, grand-canonical, 
and generalized ensembles. 



III. OVERLAPS 

We first focus on the behavior of the weights |CqP de- 
termined by the initial state and compare it with the one 
given by the canonical ensemble e^^"^^^^ JZq,^. Most of 
the results reported in this manuscript are obtained from 
calculations for supcrlattices with period two. What that 
means is that in Eq. ([T]), 

y/^' = A{-\Y. 

We will mainly focus on fillings (i) iV = L/2 (half fill- 
ing), for which observables after relaxation were seen to 
quickly approach the thermal predictions when the value 
of Ai was increased and Ap = 0, but no such thing was 
observed when Aj = Q and Ap was increased 0, and 
(ii) A'' = L/4 (quarter filling), which does not exhibit an 
approach to the thermal predictions, like the one seen at 
half filling, no matter the selected values of Ai and Ap- 
Some results for trapped systems, related to the findings 
in Ref. [21j, will be reported in the following section. 

By comparing Eqs. ([S]) and (l6|), one may naively think 
that for those states for which an approach to ther- 
mal expectation values was observed, the weights of the 
initial state in the eigenstates of the final Hamiltonian 
I Cap may approach those of the canonical ensemble 
g-E„/fcBT/2-(.g. In Fig. m we plot the values of |C„|2 
(top half in both panels) and e~^°'*''^"^/ZcE (bottom 
half in both panels) for quenches from the ground state 



of a superlattice {Aj = 4) to the homogeneous lattice 
{Ap = 0) (a) and from the ground state of the homo- 
geneous lattice {Aj = 0) to the superlattice {Ap — 4) 
(b). 

Figure [1] clearly shows that the actual values of |CqP 
differ not only quantitatively (several orders of magni- 
tude) from those of e~^°/'^^^/ZcE but also qualitatively 
different for both quenches, as the former exhibit a slower 
decay with the energy of the eigenstates. No conver- 
gence between the values of IC^p and e~^°/'^^-^/ZcE is 
observed as Aj and Ap are changed (not shown). In 
Fig. [11 we also provide information about the number 
of states, per unit area in the plot, that have a given 
weight within in each ensemble (color scale). For the 
quench from Aj = 4 to Ap = 0, one can see in Fig. \V[a.) 
(top half) that the number of states with nonzero values 
of IC'qP continuously increases as the energy increases 
and reaches a maximum around the center of the spec- 
trum, where the density of states is largest. A similar 
behavior can be seen within the canonical ensemble in 
Fig. [IKa) (bottom half). For the quenches from Aj = 
to Ap = 4, on the other hand, there are isolated islands 
with nonzero weights both in the diagonal and canon- 
ical ensembles [Fig. [l{b)]. This is because the many- 
body spectrum exhibits bands of eigenstates separated 
by gaps, which are determined by the value of Ap. Such 
a behavior can be straightforwardly understood from the 
single-particle band structure. In the periodic case, a rea- 
sonably good approximation for large systems with open 
boundary conditions, the latter exhibits two bands given 
by the expression 



e±(fc) = ±^4t2cos2(fca)-f A2, 



(19) 



where "-I-" denotes the upper band and "— " the lower 
band and k is the single particle momentum. Depending 
on which values of k are occupied in the many-body state, 
the bands seen in Fig. [Hb) form. 

Results for the same quenches as in Fig. [H but for 
quarter-filled systems, are presented in Fig. [5] The latter 
are qualitatively similar to the former in everything, ex- 
cept for the behavior of the number of eigenstates with 
nonzero values of jCcp in the quenches from Aj = 4 
to Ap = [top half in Fig. [2la)]- At quarter filling, 
when Aj — 4, the initial state imprints a modulation on 
the number of eigenstates with nonzero [Cd^ (note that 
the spectrum in the final Hamiltonian, when Ap = 0, 
has no gaps). That modulation is not present for the 
quenches at half filling [top half in Fig. [TJa)] and, as ex- 
pected (because of the continuous spectrum of the final 
Hamiltonian) , it is not present in the canonical results in 
the bottom half of Fig. [IJa). 

The fact that the weights in the diagonal and canon- 
ical (or any other) ensemble differ from each other is 
generic for integrable and nonintegrable systems [331 and, 
as such, need not preclude thermalization. After all, the 
weights with which eigenstates of the Hamiltonian con- 
tribute to the canonical and microcanonical ensembles 
also differ. The relevant quantity to compare different 
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FIG. 3: (Color online) Energy density p{E) for the quenches 
depicted in Figs. [T] and O Results are presented for the case 
Ai = 4, Af = in the left panels [(a) and (c)] and for Aj = 0, 
Ap = 4 in the right panels [(b) and (d)] and for systems at half 
filling [(a) and (b)] and at quarter filling [(c) and (d)]. p{E) 
is reported for the diagonal and canonical ensembles, and, 
when appropriate, we have fitted the results to a Gaussian 
(continuous lines in the plots). In all cases SE = 0.1. 



ensembles is the energy density p{E), which is equal to 
the sum of the weights studied in Figs. [1] and [U over a 
given energy window 6E, divided by dE. By construc- 
tion, the integral of this quantity over the full energy 
spectrum is normalized to 1. {6E needs to be selected 
in such a way that the results for the energy density are 
independent of its actual value.) The energy density de- 
pends not only on the weights but also on the density of 
states, and tells us which part of the spectrum is the one 
that contributes the most to the ensemble averages. 

In Fig. [21 we present p{E) for the quenches for which 
the weights of the diagonal and canonical ensembles were 
reported in Figs. [T] and [H As expected, the energy density 
in the canonical ensemble is very close to a Gaussian 
p{E) = (v^^£;)-ie-(^-^^)'/(25£;^) in all cases. For the 
quenches from A/ = to Ap = 4, the Gaussian is cut by 
the bands described previously. 

In diagonal ensemble, p{E) has been recently shown to 
be very well described by a Gaussian for nonintegrable 
systems and sparse (very different from Gaussian) in in- 
tegrable systems [Sa] . We find the latter to be generic for 
our quenches, as shown in Figs. [2Ic), [2Id), and, maybe 
less evident but still true. Fig. [3Kb). Surprisingly, we find 
that for quenches from Aj y^ to Ap = at half filling, 
the energy density in the diagonal ensemble approaches 
a Gaussian as the value of Aj is increased. See Fig. ^a.) 
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FIG. 4: (Color online) Scaling of the energy density in the 
diagonal (a) and canonical ensembles (b) with increasing sys- 
tem size. Results are reported for the quenches from Aj = 4 
to Af = at half filling. Continuous lines depict the result 
of the fit of each data set to a Gaussian. 



for Ai — A and Ap = 0. This highlights the special char- 
acter of this class of initial states and will be analyzed 
more quantitatively in the following sections. 

A remark is in order on the scaling of the plots shown 
in Fig. [3] with increasing system size. For the canonical 
ensemble, it is known that the width of p{E)^ relative to 
the full width of the spectrum, vanishes in the thermo- 
dynamic limit. The question is then what happens for 
the diagonal ensemble. On general grounds, for Hamil- 
tonians containing only finite-range terms, it was shown 
in Ref. i^] (supplementary materials) that the width of 
p{E), relative to the full width of the spectrum, also van- 
ishes in the thermodynamic limit. The scaling of the 
width of p(E) depends in this case on the nature of the 
quench [23]. In Fig. 21 we show a finite size scaling for 
p{E) in the diagonal (a) and canonical (b) ensembles in 
the quenches from Aj ~ 4 to Ap = 0. These results are 
consistent with the vanishing of the width of p{E), rela- 
tive to the width of the spectrum, as the system size is 
increased. 



IV. ENTROPIES 



In Fig. [5l we show the entropy per site for two differ- 
ent quenches in half-filled systems, with increasing sys- 
tem size. For both quenches, one can see that all entropy 
per site plots tend to saturate to a constant value with 
increasing L, making evident the additivity of this ob- 
servable in all ensembles. Another result that is appar- 
ent from those plots is that Sd is smaller than all other 
entropies, and it seems that it will remain that way in 
the thermodynamic limit, as noted in Ref. (36[. An im- 
portant difference between the behavior of the entropies 
for a quench from the superlattice to the homogeneous 
lattice [Fig. [SJa)] and the quench from the homogeneous 
lattice to the superlattice [Fig. ^h)] is that, in the for- 
mer, the entropy of the GGE and the grand-canonical 
ensemble are nearly identical and the entropies of the 
GCE and the canonical ensemble approach each other 
with increasing system size. In the latter quench, the 
entropies of the GGE and the grand-canonical ensemble 
differ from each other and their difference is seen to re- 
main constant as the system size is increased. S'gce and 
S'cE approach each other as the system size increases, 
but their difference is clearly larger than for the Aj ~ 2 
to Ap ^ quench. 




L L 

FIG. 5: (Color online) Entropy per site vs. L in quenches 
from the ground state of a superlattice with Aj = 2 to the 
homogeneous lattice Af = (a) and from the homogeneous 
lattice Ai — to the superlattice with Af ~ 2 (b). For 
both quenches, we show results for the diagonal Sd, canonical 
ScE, grand-canonical Sge, generalized canonical S'gce, and 
generalized grand-canonical S'gce entropies. The systems are 
at half filling N = L/2. 



In the previous section, we have shown that the en- 
ergy distribution in quenches whose initial states are the 
ground state of half-filled systems with a superlattice 
{Ai ^ 0) can be weh described by a Gaussian, typical of 
thermal states, as the value of Aj is increased. However, 
at least for the finite systems we can solve numerically, 
we showed that such a Gaussian like energy distribution 
clearly differs from that of the canonical ensemble. In 
this section, we use the entropies, including the diago- 
nal entropy Sd [3^ l3g| , as a way to quantify the scaling 
of the energy distributions in all ensembles as the sys- 
tem size is increased. In Ref. [3a], it was already shown 
that the diagonal entropy in integrable systems increases 
nearly linearly with system size, demonstrating its addi- 
tive character. 



In order to quantify the observations above for differ- 
ent quenches and fillings, in Fig. [51 we plot the scaling 
of (S'cE - Sgce)/L and {Sqe - Sgge)/L with system 
size. The left panels [Figs. [6ja) and [6]Jb)] depict the re- 
sults at half filling. Figure [HJJa) shows that for any given 
pair Ai = X -^ Ap = and A/ = — > Ap = x, where 
X = 2,4,6,8, the difference (S'cE — Sgce)/L saturates 
at greater values for the quenches starting from the ho- 
mogeneous lattice than for those starting from the su- 
perlattice (which, for the lattice sizes shown, still keep 
decreasing as the system size is increased). The differ- 
ence (S'cE — Sgge)/L, in Fig. [6jb), exhibits and even 
more remarkable behavior. It does not change with in- 
creasing system size, and it can be seen to be orders 
of magnitude smaller for the quenches starting from the 
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FIG. 6: (Color online) Difference between the entropy of the 
canonical ensemble and the GCE [(a) and (c)] and between 
the grand-canonical ensemble and the GGE [(b) and (d)] for 
quenches at half filling [(a) and (b)] and quarter filling [(c) 
and (d)]. In all panels, the results for the quenches from the 
superlattice to the homogeneous lattice are depicted using 
open symbols, while the ones from the homogeneous lattice 
to the superlattice are depicted using solid symbols. For the 
former quenches, results are reported for Ai = 2 and Af = 0, 
Ai = ^ and Af ^ 0, Ai = & and Af = 0, and ^/ = 8 and 
Af ~ 0, and, for the latter, results are reported for Ai = 
and Af ^ 2, Ai = and Af = 4, Ai = and Af = 6, and 
Ai — and Af = 8. In the legend, we use the notation "A/ 
to Af" to label the plots. 



superlattice when compared to those starting from the 
homogeneous system. The difference (S'ge ~ Sqqe)/L 
quickly approaches zero as the value Aj in the superlat- 
tice is increased. This is exactly the same behavior that 
was observed in Ref. Q for the difference between the 
expectation value of the momentum distribution func- 
tion {uk) in the grand-canonical ensemble and that of 
the time average in the time evolving state. (The latter 
can be reproduced using the GGE.) 

Hence, we can conclude that for the particular class 
of initial states in Ref. Q , where quenches starting from 
the ground state of a system with a superlattice lead 
to expectation values of n^ that approach those in ther- 
mal equilibrium as Aj was increased, the sets of states 
that contribute to the grand-canonical ensemble and the 
GGE become increasingly similar to each other. Since 
the GGE describes observables in the integrable system 
after relaxation [3, l2l[ , thermal ensembles then will also 
provide a very good estimate for those observables as Aj 
is increased. From the results in Fig. Hfb) , it is impor- 
tant to stress that the entropies in the grand-canonical 
ensemble and the GGE do not approach each other, for 
a fixed value of Aj, as the system size is increased. 

In Figs.inijc) and|ni[d), we show results for an identical 
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FIG. 7: (Color online) Difference between the entropy of 
the grand-canonical ensemble and the GGE for quenches at 
half filling (two lower curves) and quarter filling (two upper 
curves) vs. Ai. In all cases Af = 0. The dotted line depicts 
a power-law fit to the large Aj results at half filling. For each 
quench, results for two different system sizes are presented. 



set of quenches as the one in Figs.jBIa) andHl^b) but for 
systems at quarter filling. For all quenches at quarter fill- 
ing, one can see that the differences between the entropy 
in the standard ensembles and in the generalized ones is 
orders of magnitude larger than for the quenches at half 
filling. The differences between the two are maximal for 
the quenches with Aj ^ 0. The behavior with changing 
system size is, however, similar to the one in the systems 
at half filling. Hence, by comparing all panels in Fig. [51 
one can further see that there is something special about 
the quenches starting from the half-filled superlattice. 

As discussed in Refs. [a Is^l, the ground state of half- 
filled systems in a superlattice is insulating and, as the 
value of A increases, its wavefunction approaches that 
of a trivial Fock state [a product state of empty (low 
chemical potential) and occupied (high chemical poten- 
tial) sites] . In Ref. [^ , it was shown that the one-particle 
correlation length ^ decays as a power law ^/a ^ Xj ^ Ajt 
for large values of Ajt [A/t > 4). In Fig. [3 we show how 
(S'ge — Sgge)/L decreases as Aj increases. Here again, 
we find a power- law decay for large values of A/, where 
(S'ge ~ Sgge)/L ^ 1/^/- This large exponent explains 
the fast reduction of (S'ge ~ Sqqe)/ L seen in Fig. IHlwhen 
A J was increased. 

Since for large values of Aj the initial states are nearly 
uncorrelated states {^ — ?> 0), their overlaps with the 
eigenstates of the final Hamiltonian can be understood 
to be random and constrained only by energy conserva- 
tion. This helps in understanding the origin of Gaus- 
sian energy distribution observed in the previous section 
and the closeness of the generalized ensemble entropies 
to those of the standard ensembles as Aj is increased. In 
Fig. ini we also present results for the quenches at quar- 
ter filling, where (Sge — Sgge)/^ is seen to saturate to 
a finite value when Aj is increased. For both fillings, and 
the system sizes depicted in that figure, finite-size effects 
can be seen to be negligible. 

Confirmation of the conclusions above can be obtained 
if one realizes that a similar argument applies to the sys- 




FIG. 8: (Color online) (a) Density in the center of the trap 
as a function of the excitation energy per particle e, which is 
changed by increasing Vi in a system with 50 lattice sites and 
10 particles, (b) Difference between the entropy of the grand- 
canonical ensemble and the GGE vs. L for systems with 
different excitation energy per particle. (Inset): Integrated 
relative difference between rik in the diagonal and canonical 
ensembles (see text) vs. the excitation energy per particle. 
Results are presented for different system sizes, L = 25, 30, 
...,45 fH. 



terns discussed in Ref. [2l[. There, the initial state was se- 
lected to be the ground state of a trapped system, where 
[Eq. (HD] 

F/'^* = V{j - L/2f 

is a harmonic trapping potential and the evolution was 
followed after the trap potential V was turned off, i.e., 
V7 7^ and Vp = 0. For a fixed number of particles, 
as Vi increases, a Mott insulator (Fock state for hard- 
core bosons) with density n = I forms in the center of 
the trap. When initial states containing such Mott in- 
sulating domains were used for the time evolution, the 
difference between the momentum distribution function 
in the diagonal ensemble and standard ensembles of sta- 
tistical mechanics was seen to decrease |2ll |. 

In the inset in Fig. [SJ we show the results obtained in 
Ref. [2l|] for the integrated difference between the predic- 
tions of the diagonal and canonical ensembles for n^ 
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as a function of the excitation energy per particle 

Ej — Eq 

where Eg is the ground-state energy of the final (homo- 
geneous) Hamiltonian. The excitation energy increases 
by increasing Vj, while keeping L and N constant |2l| . 

The density in the center of the trap (in the initial 
state) versus the excitation energy is plotted in Fig. M,^) ■ 
There, one can see that (Anfe)cE (in the inset) is smallest 
and keeps decreasing when the density in the center of 
the trap approaches or becomes equal to 1, i.e., when an 
increasingly large portion of the system comes close or 
becomes a Fock state. 



The scaling of the difference between the entropies 
in grand-canonical ensemble and the GGE is shown in 
Fig. m^b) for different excitation energies. Similarly to 
the results for the superlattice systems, that difference is 
seen to be smallest (and decreasing with increasing sys- 
tem size in this case) for the systems whose initial states 
are closest to Fock states. Hence, once again, a spe- 
cial class of initial states is seen to produce increasingly 
"thermal-like" observables and generalized ensembles. 



V. CONSERVED QUANTITIES 

Conserved quantities play a fundamental role in the 
dynamics and description after relaxation of integrable 
systems. The latter follows from the evidence that gen- 
eralized ensembles are able to describe observables after 
equilibration while standard statistical ensembles are, in 
general, not [Ml [iMl, HH-il. Hence, a distinctive 
behavior is expected of the distribution of the conserved 
quantities for those initial states for which observables 
after relaxation approach thermal values. In this section, 
we study the behavior of the conserved quantities in those 
and other cases analyzed in the previous sections. 

As explained in Sec. HIl the expectation values of the 
conserved quantities in hard-core boson systems can be 
straightforwardly computed because they are the occu- 
pation of the eigenstates of the noninteracting fermionic 
Hamiltonian (there are L of those) to which hard-core 
bosons can be mapped. As such, they can be calcu- 
lated within the GGE (identical to those of the initial 
state and the diagonal ensemble by construction) and in 
the grand-canonical ensemble, for very large lattices. In 
the grand-canonical ensemble, the occupation of the con- 
served quantities is dictated by the Fermi distribution 
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where e„ are the single-particle eigenenergies. 

In Figs. El^a) and|9l^b), we depict the conserved quan- 
tities in the GGE (initial state) and the grand-canonical 
ensemble for quenches at half filling from the ground state 
in a superlattice to the homogeneous lattice [Fig. EJa)] 
and from the ground state of the homogeneous lattice 
to the superlattice [Figs. Eljb)]. (The conserved quan- 
tities are ordered from the highest to the lowest occu- 
pied in the initial state.) A clear contrast can be seen 
between those two panels. In Fig. [9l[a) the results for 
the GGE and grand-canonical ensemble are almost in- 
distinguishable from each other while in Fig. ^h) they 
differ from each other markedly. This behavior docs not 
change with increasing system size as, in the same fig- 
ure, depicted as continuous lines, we also report results 
for lattices 10 times larger than those used for the calcu- 
lations depicted as symbols. Qualitatively, these results 
are similar to those obtained in Ref. [21| for the trapped 
systems analyzed in the previous section. 
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FIG. 9: (Color online) Expectation value of the conserved 

quantities in quenches from Aj j^ to Af = [(a) and (c)] 

and from Ai — to Af 7^ [(b) and (d)] for systems at half 

filling [(a) and (b)] and systems at quarter filling [(c) and (d)]. 

The conserved quantities are ordered from the highest to the 

lowest occupied in the initial state. Solid symbols depict the 

results of the GGE (conserved quantities in the initial state) 

and open symbols depict the results of the grand-canonical 

ensemble. The results denoted by symbols (lines) correspond 

to systems with 38 (380) sites in the half-filled case [(a) and 

(b)] and with 48 (480) sites in the quarter-filled case [(c) and 

(d)]. Note that for the smallest system sizes depicted here (the 

largest analyzed in the previous sections) finite-size effects for 

the conserved quantities already are negligible. They exhibit 

an almost perfect overlap with the results in systems 10 times 

larger. Results are reported for quenches between Aj — 2 and 

Af ^0, Ai = A and Af = 0, Ai = 6 and Af = 0, and Ai = 8 

and Af = and between Aj — and Af = 2, Aj = and 

Af = 4:, Ai ^0 and Af = 6, and Ai = and Af = 8. In the 

legend, we use the notation "Ai to Af" to label the plots. 



Insights into the contrast between the results in 
Fig. M,&) and Fig. MJ^) can be gained if one notices that 
the distribution of conserved quantities for the quenches 
from and to the superlattice are smooth and identical 
when Aj in the former is equal to Ap in the latter. 
This immediately helps one understand why the grand- 
canonical ensemble prediction of the conserved quanti- 
ties can match that of the quenches from the superlat- 
tice {Aj 7^ 0) to the homogeneous lattice (Ap = 0) but 
not that of the quenches from the homogeneous lattice 
{Aj = 0) to the superlattice {Ap 7^ 0). In the former, the 
final system exhibits no gaps (Ap = 0) and the conserved 
quantities [the Fermi distribution, see Eq. (pO)) ] can be a 
smooth function of n at finite temperatures, while in the 
latter the system is gapped {Ap ^ 0) and, hence, a dis- 
continuity must occur in the Fermi distribution at the 
gap position [as seen in Fig. Eljb) for n/L = 0.5]. 

The results reported in Figs. [5Jc) and Eljd) for sys- 
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FIG. 10: (Color online) Integrated differences between the 
conserved quantities in the GGE (initial state) and in the 
grand-canonical ensemble (see text) vs. L for half-filled sys- 
tems in quenches from Ai ^Q to Af ~ (a) and from Ai — 
to Af 7^ (b). Results are reported for quenches between 
Ai ^ 2 and Af = 0, Ai = 4 and Af = 0, Ai = 6 and 
Af = 0, and Ai — 8 and Af = in (a); and, between Ai — 
and Af = 2, Ai = and Ap = 4:, Ai = and Ap = 6, and 
Ai = and Af = 8, in (b). In the legend, we use the notation 
"Ai to Ap" to label the plots. 



tcms at quarter filling, are qualitatively similar to those 
in Fig.[9l^b). For all quenches, the conserved quantities in 
the initial state differ substantially from those predicted 
by the grand-canonical ensemble. The differences can 
be noted to be particularly large if one realizes that, for 
many conserved quantities, the initial state has a zero ex- 
pectation value while the grand-canonical ensemble pre- 
dicts nonzero, and large, values. This helps in under- 
standing the large differences seen in the previous section 
between the entropies in the generalized and standard en- 
sembles for the quenches at quarter filling. Once again, 
the behavior of the conserved quantities in thermal equi- 
librium is dictated by the Fermi distribution and can be 
understood given the gaplcss or gapped nature of the 
spectrum of the final system. 

Figure [10] depicts how the difference between the con- 
served quantities in the initial state and the grand- 
canonical ensemble, given by the integrated relative dif- 
ference 
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behaves as the system size increases. The results pre- 
sented, for half-filled systems in quenches from a super- 
lattice potential {Aj 7^ and Ap = 0) in Fig. [TOla) 
and to a superlattice potential {Aj ~ and Ap 7^ 0) in 
Fig. llOf b). show more quantitatively that the results in 
Fig. [S] do not change with increasing system size. 

Figure [To] also makes evident that there is a big quanti- 
tative difference between A/ in the systems whose initial 
state is the ground state in the superlattice [Fig. (TUT a)] 
and those whose initial state is the ground state of the 
homogeneous lattice [Fig. (TOT b)]. This is better seen in 
Fig. [Til where A/ is plotted versus Aj for the former 
case and versus Ap for the latter. In both cases, we 
find power-law decays, which are ^ 1/j4j when the ini- 
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FIG. 11: (Color online) Integrated differences between the 
conserved quantities in the GGE and in the grand-canonical 
ensemble. Results are reported for quenches where (i) Aj j^ 
and Af ~ (two bottom curves) vs. Ai and (ii) ^4/ = 
and Af 7^ (two upper curves) vs. Af and for two different 
system sizes. The dotted lines depict a power-law fits to the 
large Ai, Af results. 



tial state was created for Aj y^ and ^ l/Ap when the 
final Hamiltonian has Ap 7^ 0. It is important to stress 
that while increasing Ai does not qualitatively change 
the time dynamics of observables of interest, increasing 
Ap does @. In the latter case the damping (relaxation) 
of the observables is inhibited @, so the assumption 
that observables relax to time independent values breaks 
down. 

Finally, from Figs. [TU] and [TTl we should emphasize 
once again that, complementary to the behavior seen for 
(•Sge ~ Sqq-e)/L in the previous section, the scaling of 
A/ versus L is similar for both types of quenches, namely, 
any finite value of Aj (if Ap = 0) or Ap (if Aj = 0) leads 
to a finite A/ in the thermodynamic limit. The difference 
between those quenches resides in the actual values of A/ 
and their behavior with changing Aj or Ap. 

On the basis of those results we can now understand 
that, for the classes of initial states in Refs. [^, [21| for 
which observables after relaxation approached thermal 
values, the control parameter used tuned the distribu- 
tion of conserved quantities to approach thermal values 
(resulting in generalized ensembles that, for those states, 
approach thermal ensembles). This behavior, however, 
should not be confused with thermalization as under- 
stood for nonintegrable systems. For the latter, the dif- 
ference between observables after relaxation and the pre- 
dictions of statistical mechanics ensembles is expected to 
vanish in the thermodynamic limit, while, for the special 
classes of initial states that we have studied here for in- 
tegrable systems, such a difference remains finite in the 
thermodynamic limit for any selected (finite) value of the 
control parameter. 

To conclude, there is an important distinction to be 
made about the generalized ensembles when compared 
with standard ensembles of statistical mechanics. In the 
latter, the conserved quantities (energy, momentum, an- 
gular momentum, etc.) are additive and their number is 



^ 1. In the generalized ensembles, the conserved quan- 
tities are, strictly speaking, not additive and their num- 
ber, in the integrable systems considered here, is ^ L. 
The fact that they are not additive can be immediately 
seen in Fig. [9l where, after making the system size 10 
times larger, the value of the conserved quantities does 
not change. Instead, their number increased by a fac- 
tor of 10 (that is the reason for plotting the conserved 
quantities as functions of n/L). 
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FIG. 12: (Color online) Lagrange multipliers in quenches from 
Ai ^ Q to Af ~ for systems at half filling (a) and systems 
at quarter filling (b). (The results for quenches from Aj = 
to Af 7^ are the identical.) Symbols (lines) correspond to 
systems with 38 (380) sites in the half-filled systems [(a) and 
(b)] and to systems with 48 (480) sites in the quarter-filled 
systems. Once again, note that size effects for the Lagrange 
multipliers are negligible. Results are reported for quenches 
between Ai = 2 and Af = 0, Ai = A and Af = 0, Ai = 6 
and Af = 0, and Ai = 8 and Af = 0. In the legend, we use 
the notation "Ai to Af" to label the plots. 

In Fig. [T21 we show the values of the Lagrange multi- 
pliers for the same quenches and system sizes depicted in 
Fig. ini As expected from the expression for the Lagrange 
multipliers [Eq. ((T^ ]. they are a smooth function of the 
values of the conserved quantities (and exhibit negligible 
finite-size effects in Fig. [T2)) . One can then think of the 
conserved quantities, considered here to build the gener- 
alized ensembles, as additive in a coarse-grained sense. 
This follows if one realizes that by, increasing the system 
size, the Lagrange multipliers in a coarse-grained region 
do not change their values (Fig. [T^ . but the sum of the 
expectation values of the conserved quantities in that 
region (Fig. |9]) grows proportionally to the increase of 
system size. Hence, effectively, the conserved quantities 
behave as additive. A discussion on the role of additivity 
of the conserved quantities in generalized ensembles can 
be found in Ref. uM- 



VI. SUMMARY 

We have studied the dependence on the initial state 
of the description of integrable systems after relaxation 
following a sudden quench. In general, integrable sys- 
tems are not expected to thcrmalize. Hence, we have 
focused on understanding special classes of initial states 
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that lead to values of observables after relaxation that 
approach those in thermal equilibrium, when a control 
parameter is changed. One of our main findings is that, 
even for such initial states, thermalization does not oc- 
cur as in nonintegrable systems. In the latter, the differ- 
ence between the thermal expectation value of an observ- 
able and those after relaxation is expected to vanish in 
the thermodynamic limit. In the intcgrable systems dis- 
cussed here, no matter the initial state selected (which 
is an eigenstate of another integrable system where the 
control parameter is one of the parameters of the initial 
Hamiltonian) , the distribution of conserved quantities in 
the thermal ensembles differs from (but can be arbitrar- 
ily close to) that of the diagonal ensemble (or the GGE) , 
and the difference does not vanish with increasing sys- 
tem size. Since the values of the conserved quantities 
constrain the outcome of the relaxation dynamics, the 
observables after relaxation do not reach thermal values 
in the thermodynamic limit. 

Another of our main findings is that what the con- 
trol parameter is doing in those special classes of initial 
states is tuning the distribution of conserved quantities to 



approach thermal values. As a result, the initial states 
exhibit energy densities that are increasingly Gaussian 
like and entropies of their associated generalized ensem- 
bles that approach those of standard ensembles. Sim- 
ilarly to the behavior seen for the conserved quantities, 
the difference between the entropy per site in the general- 
ized and standard ensembles remains nonzero in the ther- 
modynamic limit. It can, however, be made arbitrarily 
small by changing the control parameter. Interestingly, 
for the model considered here, the special initial states 
were found to be insulating ground states that approach 
products of single site wavefunctions. 
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